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SUMMARY 

Large-displacement elastic and elastic-plastic, finite-element stress-strain analyses of an oxygen- 
free high-conductivity (OFHC) copper plate specimen were performed using an updated Lagrangian for- 
mulation. The plate specimen is intended for low-cost experiments that emulate the most important 
thermomechanical loading and failure modes of a more complex rocket nozzle. The plate, which is loaded 
in bending at 593 °C, contains a centrally located and internally pressurized channel. The cyclic crack 
initiation lives were estimated using the results from the analyses and isothermal strain-controlled low- 
cycle fatigue data for OFHC copper. A comparison of the predicted and experimental cyclic lives showed 
that an elastic analysis predicts a longer cyclic life than that observed in experiments by a factor greater 
than 4. The results from elastic-plastic analysis for the plate bend specimen, however, predicted a cyclic 
life in close agreement with experiment, thus justifying the need for the more rigorous stress-strain 
analysis. 


INTRODUCTION 

One of the important life-limiting failure modes encountered by the rocket nozzle liner of the 
Space Shuttle Main Engine (SSME) is thermal/mechanical ratchetting of the thin (0.0635-cm) wall of 
copper material between the fuel/coolant channels and the hot gas wall. The ratchetting of material is 
toward the inside hot gas path of the rocket nozzle. When viewed in transverse (to the gas flow) cross 
section, the initially rectangular fuel/coolant channels tend to develop a “doghouse” shape with the roof 
pointed toward the centerline of the nozzle (see fig. 1 from ref. 1). This deformation mode feeds on itself 
and causes the thin wall to become ever thinner. The thinned area may rupture eventually and dump valu- 
able hydrogen fuel into the nozzle where it is burned without a significant propulsive effect. The rate of 
thinning (i.e.. ratchetting) has been shown analytically by Arya and Arnold (ref. 2) to be temperature, 
pressure, and cycle dependent. New geometric designs or new materials are required to either eliminate or 
drastically reduce this failure mode. 

Full-scale nozzle liners are expensive to build and test for the purpose of investigating new con- 
figurations and construction materials. Thus, subscale liners that retain the important features of a 
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full-scale liner have been designed, built, analyzed, and tested at the NASA Lewis Research Center 
(refs. 3 and 4). These too are expensive to build and test. Furthermore, the turnaround time from idea 
(new material, new channel configuration, fabrication, etc.) to experimental verification is unacceptably 
long (typically about 2 yr). Thus, if a simpler, less expensive benchmark configuration can successfully 
emulate the basic structural loading and failure mode of full-size hardware, tremendous savings can be 
achieved. 

A relatively simple (ground) plate bend specimen developed at Lewis (LJ. Westfall, 1992, NASA 
Lewis Research Center, Cleveland, OH, private communication ) is able to capture the essence of the load- 
ings and failure modes found in high-heat-flux, high-thermal-conductivity rocket nozzle liners such as the 
Space Shuttle Main Engine (SSME). The bend specimen is a flat plate with an internally pressurized 
single longitudinal channel that is located along the center of the plate (fig. 2). The appendix gives the 
manufacturing details of this specimen. 

Large-displacement tests of the plate specimen were performed isothermally at a high temperature 
in four-point loading to simulate the large cyclic thermal strains experienced by the full-scale nozzle liner. 
Bending strain gradients in the plate also simulate thermal strain gradients in the nozzle liner. That this 
bend specimen emulates the failure mode observed in nozzle liners is seen in figure 3 where the doghouse 
effect is evident: the deformed configuration of the bend specimen observed in tests performed at Lewis. 
Therefore, we believe that this specimen meets the requirements of a less expensive benchmark configu- 
ration that can be used to analytically assess and experimentally verify new materials and design concepts 
for reducing ratchetting. 

This paper analytically investigates the loading and failure modes in the surrogate bend specimen. 
Because the ratio of the maximum deflection of the plate to its thickness is large (greater than 0.75), 
small-displacement analyses for the problem produced inaccurate results. Large-displacement analyses 
are therefore required. Although the temperature (593 °C) of the copper plate is in its creep range, the 
testing frequencies are high enough (6.5 cycles/min) to justify finite-element analyses in which the time 
(creep) effects are neglected. The investigations reported herein were conducted by performing large- 
displacement finite-element analyses in conjunction with elastic and elastic-plastic constitutive models. 
The analytical results are displayed and discussed. The future direction of this research work is also 
discussed. 


FINITE-ELEMENT ANALYSES 
Finite-Element Model 

The geometry and dimensions of the specimen and channel are given in figure 2. Solid eight- 
noded brick elements with an assumed strain formulation were utilized to model the pressurized 
(6.2-MPa) center channel of the specimen. Because of the geometric and loading symmetry, only one- 
quarter of the specimen was modeled. A sensitivity study of the effect of mesh density on the results 
showed that a large number of elements was required for accurately capturing the high bending strain gra- 
dients. The presence of internal pressure in the channel, in addition to the bending gradients, necessitated 
the use of an even finer mesh in the channel region of the specimen. This finite-element model is shown 
in figure 4 and consists of 3330 elements and 4463 nodes. 


Approximate Shape of Bent Specimen (Plate) 

Only the maximum deflection occurring at the center of the plate was experimentally measured; 
hence, to determine the deflections at all other points of the plate (and its bent shape at maximum deflec- 
tion), an approximate elastic analysis was performed by modeling the whole plate as a solid thick plate 


2 


excluding the channel under uniform load transverse to the plate. The load was analytically increased 
until the calculated deflection at the center was equal to the value measured in experiments. The deflec- 
tions at all other points (nodes) of the plate were also recorded at this time. These displacements were 
later used to perform the elastic and elastic-plastic finite-element analyses by treating these latter analyses 
as displacement-controlled boundary value problems. All finite-element analyses were performed by 
using the finite-element program MARC (ref. 5). The undeformed and (calculated) deformed shapes of 
the solid thick plates subject to these finite-element analyses are shown in figure 5. The displacements 
shown are exaggerated by a factor of 5. 


Elastic Analysis 

A one-fourth segment of the pressurized and deformed (bent) single-channel plate is shown in fig- 
ure 6. The undeformed configuration is presented for comparison. Displacements are exaggerated by a 
factor of 5 again to clearly distinguish the two configurations. An elastic stress-strain structural analysis 
was performed using a finite-element discretization of this segment as described in the section Finite- 
Element Model. Because the problem involves large displacements of the segment in the XZ-plane, a 
large-displacement analysis using an updated Lagrangian formulation was adopted. The channel region of 
the specimen is of greatest interest; thus, only those figures that exhibit the stress and strain distributions 
in this region are included. 

The longitudinal (X-direction) stress in the vicinity of the channel is shown in figure 7. A typical 
bending stress distribution is observed with tension on the upper surface, compression on the lower sur- 
face, and a neutral axis of zero stress in the middle. The maximum elastically calculated longitudinal 
stress in the channel is about 766 MPa. The stresses in the Y- and Z-directions are shown in figures 8 
and 9, respectively. The Y- and Z-stresses have maximum magnitudes of about 350 and 260 MPa, 
respectively. 

The strain distribution in the channel in the X-, Y-, and Z-directions is shown in figures 10 to 12, 
respectively. The maximum magnitudes of the strains occur in the X-direction (approx equal to 0.7 per- 
cent). Also seen in figure 10 is the X-strain, which is tensile in the top layer of the channel and compres- 
sive in the bottom layer. The middle of the channel has almost zero X-strain. 


Elastic-Plastic Analysis 

Elastically calculated stress values far exceed the 0.2-percent offset yield strength of OFHC cop- 
per (30.0 MPa) at 593 °C. Therefore, an elastic analysis of the problem is inappropriate and may lead to 
erroneous results. Obviously, inelastic constitutive behavior must be taken into consideration. Keeping 
this in mind, a large-displacement elastic and perfectly plastic finite-element analysis of the segment was 
performed. As in the case of the elastic analysis, the large-displacement elastic-plastic analysis is also 
based on an updated Lagrangian formulation. Some important results from elastic-plastic analysis are pre- 
sented here. 

The X-stress distribution in the channel obtained from an elastic-plastic analysis is shown in fig- 
ure 13. The maximum X-stress of 46.4 MPa contrasts with 766 MPa obtained from an elastic analysis. 

The Y- and Z-stress distributions in the channel are shown in figures 14 and 15. respectively. The maxi- 
mum Y- and Z-stress magnitudes calculated from an elastic-plastic analysis are 34.2 and 14.9 MPa, 
respectively. The stress values obtained from an elastic-plastic analysis are significantly more realistic 
than those obtained from the elastic analysis. 

The total X-, Y-, and Z-strain distributions in the channel are shown in figures 16 to 18. As 
expected, the X-strain is tensile in the top layer and compressive in the bottom layer. Similarly, because 
of Poisson’s effect, the Z-strain is compressive in the top layer and tensile in the bottom layer (fig. 18). 
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Figure 1 9 shows only the plastic component of the X-strain in the channel. A comparison of this figure 
with figure 16 reveals that the deformation of the channel is largely plastic. 


Life Analysis 

The cyclic crack initiation life of the bend specimen was estimated by using the results from large 
elastic and large elastic-plastic structural analyses. Following the procedure of Manson and Halford 
(ref. 6), the maximum total strain ranges at the critical locations of the bend specimen were estimated by 
using the following expressions: 

effective stress: 


and 




1/2 


effective strain: 







1/2 


The values of maximum total strain range for the bend specimen obtained from the elastic and elastic- 
plastic analyses are listed in table I. 


TABLE I.— EXPERIMENTAL AND PREDICTED 
CYCLIC LIVES OF PLATE BEND SPECIMEN 


Analysis 

Maximum total effective 

Life, 


strain range. 

cycles 


percent 


Large elastic 

i.i 

500 

Large elastic-plastic 

1.7 

190 

Experiment 

— 

109 


Conway, Stentz, and Berling (ref. 7) reported isothermal, strain-controlled, low-cycle fatigue test 
results for annealed, oxygen-free, high-conductivity copper at 538 °C in a nonoxidizing argon atmo- 
sphere. These results are shown in figure 20 as total effective strain range versus cycles to crack initiation. 
The data were generated at a cyclic strain rate of 0.002 s _ I . Because this strain rate is about the same 
(half) as that experienced by the bend specimen, the data of reference 7 at 538 °C can reasonably be 
applied to the assessment of the cyclic life for the bend specimen. Outputs from the elastic and elastic- 
plastic analyses give the total maximum effective strain ranges at the critical locations of the specimen as 
1.1 and 1.7 percent, respectively (see table I). From figure 20, these strain ranges are shown to correspond 
to cyclic lives of 500 and 190. respectively. Note that the maximum total effective strain ranges were 
determined from the results of elastic and elastic-plastic finite-element analyses performed at 593 °C as 
were the tests for the OFHC copper specimens. The experimental data reported in reference 7, however. 
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were obtained at 538 °C. Because no experimental data from OFHC copper are available at 593 °C, the 
538 °C data were used to assess the cyclic life for the specimen. The justification for using the experimen- 
tal data at 538 °C to predict the cyclic lives from the results of analyses performed at 593 °C lies in the 
fact that the tensile ductility of OFHC copper is insensitive to temperature (variations measured at Lewis 
range from an 81- to 89-percent reduction of area). The low-cycle fatigue resistance is, therefore, also 
expected (ref. 8) to be insensitive to temperature over the range of interest. 

Westfall (1992, NASA Lewis, Cleveland, OH, private communication) experimentally determined 
the cyclic life of the bend specimen to be 109 at 593 °C. The elastic analysis predicts a longer cyclic life 
than that observed in experiments by more than a factor of 4. Although the elastic-plastic analysis also 
predicts a longer cyclic life for the bend specimen, the agreement between the prediction and experiment 
is much closer. Because the life predictions are based on a pure fatigue failure mode (exclusive of 
ratchetting), it is not surprising that the predictions exceeded the experimentally observed cyclic lives. 


SUMMARY 

Large-displacement elastic and elastic-plastic finite-element analyses were presented for a grooved 
plate bend specimen. The analyses incorporate an updated Lagrangian formulation. Using the results from 
these elastic and elastic-plastic finite-element analyses, the cyclic lives for crack initiation are assessed for 
the bend specimen. The elastic analysis predicted a longer cyclic life for the specimen than that observed 
in experiments by a factor of more than 4. The elastic-plastic analysis also predicted a longer cyclic life 
than that observed, but the prediction was in closer agreement with the experimentally determined cyclic life. 


FUTURE WORK 

Large-deformation elastic analysis predicts a longer cyclic life for a grooved specimen than that 
observed in experiments. In contrast, elastic-plastic analysis predicts cyclic life much closer to the experi- 
mental result. Time (creep) effects were assumed insignificant and were excluded from the analyses. It is 
intended that the life prediction methodology described in this report be applicable to the SSME rocket 
nozzle liner. The period of a laboratory specimen test cycle is about 9 sec. The period of thermomechani- 
cal loading cycle experienced by the SSME rocket nozzle liner is more than 480 sec. An accurate assess- 
ment of the structural response of a bend specimen that emulates the realistic thermomechanical loading 
cycle experienced by an SSME nozzle liner will require time-dependent inelastic analysis. In this regard, 
viscoplastic models provide more accurate descriptions and assessments of inelastic response than those 
of conventional elastic-plastic (creep) models. It is therefore proposed that a 480-sec loading cycle- 
imposed viscoplastic finite-element analysis be performed for the bend specimen for future cyclic life 
assessment. Furthermore, to assess the creep-fatigue failure mode life for loading cycles experienced by 
the SSME rocket nozzle liners, laboratory material tests with durations of the order of 480 sec/cycle 
would be essential. 
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APPENDIX— FABRICATION OF COPPER PLATE SPECIMENS 


Oxygen-free high-conductivity (OFHC) copper sheets, 0.3175 and 0.0889 cm thick, were cut into 
5.08- by 15.24-cm sections to match the 5.08-cm-wide channel die used for specimen fabrication. The 
thicker sheets were machined to produce a center chamber, 0.3 175 by 5.08 cm, with a thin gas transport 
passage (shown in fig. 2). A 10.16-cm-long, 0.3175-cm-diameter stainless steel tube was inserted 5.08 cm 
into the gas transport passage to produce a gas pathway to the center chamber. 

Prior to assembly operations, all surfaces of the copper sheets were mechanically cleaned with fine 
sandpaper to remove any copper oxide. Nickel-palladium (Ni-35.32Pd-9.63Cr-4.63Fe2.69B) brazing foil, 

0.00381 cm thick, was used for pressure brazing the two sheets together. This brazing alloy melts at a 
temperature between 945 and 996 °C. Fifteen-cm-long sheets of the foil were cut and machined to pro- 
duce a hole that matched the center chamber of the thick copper sheet. The brazing foil was secured to 
both surfaces of the thick copper sheet by a liquid brazing flux which acted as an adhesive. A debonding 
agent or AI 2 O 3 wool was placed in the center chamber to prevent the flow of brazing material into the 
chamber. Thin copper sheets were placed on top of the brazing foil and secured with tape or copper rivets 
at the four comers. The three-sheet assembly was covered with a release agent and placed into a vacuum 
hot press containing a 5.08-cm-wide channel die. The assembly was vacuum hot pressed (pressure 
brazed) by heating the die containing the sample to 982 °C under a pressure of 3.4475 MPa in a vacuum 
of 5X10-5 torn for a hold time of 1 hr. 
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Figure 1 . — Thrust chamber coolant channels before and 
after deformation (ref. 1). (a) Before, (b) After. 




Figure 3. — Deformed channel of plate specimen showing 
"doghouse" effect. 
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Deformed 



Figure 4. — Finite-element model of plate bend specimen (3330 
elements; 4463 nodes). 



Figure 5. — Undeformed and deformed shapes of thick plate, 
(Deformations are magnified by a factor of 5.) 


Deformed 




Undeformed 


Figure 6. — Deformed and undeformed shapes of segment. 
(Deformations are magnified by a factor of 5.) 
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Figure 7. — X-stress in channel wall; elastic analysis. 
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Figure 8. — Y-stress in channel wall; elastic analysis. 
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Figure 9. — Z-stress in channel wall; elastic analysis. 
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Figure 10.— X-strain in channel wall; elastic analysis. 
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Figure 1 1 . — Y-strain in channel wall; elastic analysis. 
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Figure 1 2. — Z-strain in channel wall; elastic analysis. 
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Figure 13. — X-stress in channel wall; elastic-plastic analysis. 



Figure 14. — Y-stress in channel wall; elastic-plastic analysis 
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Figure 15 — Z-stress in channel wall; elastic-plastic analysis. 





Figure 1 6 — Total X-strain in channel wall; elastic-plastic analysis. 
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Figure 18. — Total Z-strain in channel wall; elastic-plastic analysis. 
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Figure 19. — Plastic X-strain in channel wall; elastic-plastic analysis. 
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Cycles to crack initiation 

Figure 20. — Low-cycle fatigue resistance of oxygen-free, high- 
conductivity (OFHC) copper at 538 °C in argon at cyclic strain 
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